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Abstract 


We present high-resolution infrared spectroscopic measurements of the ecliptic comet 2P/Encke, 
observed on 4-6 Nov. 2003 during its close approach to the Earth, using the Near Infrared 
Echelle Spectrograph on the Keck II telescope. We present flux-calibrated spectra, production 
rates, and mixing ratios for H 2 O, CH 3 OH, HCN, H 2 CO, C 2 H 2 , C 2 H 6 , CH 4 and CO. Comet 
2P/Encke is a dynamical end-member among comets because of its short period of 3.3 years. 
Relative to “organics-normal” comets, we determined that 2P/Encke is depleted in HCN, H 2 CO, 
C 2 H 2 , C 2 H 6 , CH 4 and CO, but it is enriched in CH 3 OH. We compared mixing ratios of these 
organic species measured on separate dates, and we see no evidence of macroscopic chemical 
heterogeneity in the nucleus of 2P/Encke, however, this conclusion is limited by sparse temporal 
sampling. The depleted abundances of most measured species suggest that 2P/Encke may have 
formed closer to the young Sun, before its insertion to the Kuiper belt, compared with “organics- 
normal” comets - as was previously suggested for other depleted comets (e.g. C/1999 S4 
(LINEAR)). We measured very low rotational temperatures of 20 - 30 K for H 2 0, CH 3 OH and 
HCN in the near nucleus region of 2P/Encke, which correlate with one of the lowest cometary 
gas production rates (~ 2.6 x 10 27 molecules s' 1 ) measured thus far in the infrared. This suggests 
that we are seeing the effects of more efficient radiative cooling, insufficient collisional 
excitation, and/or inefficient heating by fast H-atoms (and icy grains) in the observed region of 
the coma. Its extremely short orbital period, very low gas production rate, and classification as an 
ecliptic comet, make 2P/Encke an important addition to our growing database, and contribute 
significantly to the establishment of a chemical taxonomy of comets. 
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1. INTRODUCTION 


The establishment of a taxonomic classification for comets based on their chemical composition 
is essential to understanding the formation and evolution of our Solar System (for a recent 
review, see Mumma & Chamley 2011). Comets are relatively unaltered remnants from the 
creation of the Solar System 4.6 billion years ago, and are currently found in two major 
dynamical reservoirs: the Kuiper belt and the Oort cloud (Gladman 2005). 

The dynamical classification of a comet is based on its orbital Tisserand invariant with respect to 
Jupiter (Tj). For Jupiter-family comets: 2 < Tj < 3; for Chiron-type comets: Tj > 3 & a > aj; for 
Encke-type comets: Tj > 3 & a < aj; and for nearly isotropic (i.e., Oort cloud) comets: -2 < Tj < 2, 
where “a” is the semi-major axis of a comet (Levison 1996). Identifying the long-term dynamical 
reservoir of a comet does not, however, identify its formative region. The “ Nice ” model predicts 
significant dynamical dispersion of icy bodies in the outer proto-planetary disk (Tsiganis et al. 
2005), and while it suggests that comets currently found in the Kuiper belt formed beyond 5 AU 
from the Sun, the newer “Grand Tack” model (Walsh et al. 201 1) points to the possible origin of 
many Kuiper belt and Oort cloud comets in the terrestrial planets region. Another new model 
(Levison 2010) questions whether all Oort cloud comets indeed formed in our proto-planetary 
disk, and suggests that many may have been captured after ejection from the proto-planetary 
disks of other stars in the Sun’s birth cluster. Clearly, a different metric is needed to clarify a 
comet’s formative region. 

Current models for the protoplanetary disk suggest the presence of strong radial gradients in 
chemistry and temperature, and direct sampling of cometary material confirms the presence of 
pre-cometary grains that formed over extremely diverse conditions, ranging from the interstellar 
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medium to the near-Sun inner planetary region (Brownlee et al. 2006). These models also predict 
compositional and isotopic diversity amongst the ices formed in different regions, and thus in 
comet nuclei that accumulated from them. The recent discovery of ocean-like water in the 
ecliptic comet 103P/Hartley 2 (Hartogh et al. 2011) emphasizes this diversity, re-opens the 
question of the role of comets in delivering Earth’s water, and makes the need for a robust 
compositional taxonomy even more compelling. 

Optical surveys of product species (such as OH, CN, NH, C2, C3, etc.) revealed that the fraction 
of comets depleted in C2 and C3 radicals (relative to CN) is greater among Jupiter family comets 
than among Oort cloud comets (A’Heam et al. 1995, Fink 2009, Langland-Shula & Smith 20 1 1 ,a 
Cochran et al. 2012). However, it is often problematic to uniquely interpret the origin of product 
species. They can originate from the photo-dissociation of primary (a.k.a. parent) volatiles (i.e., 
species native to the cometary nucleus and released under the influence of solar radiation), from 
refractory grains, or even from chemical processes in the coma. Thus, there are multiple possible 
(and unknown) pathways and parents for each observed radical. Also, comets that are termed 
‘typical’ in product species (C2, C3, CN) are sometimes not “organics-normal” in their primary 
volatile composition (8P/Tuttle, Bonev et al. 2008; lOP/Tempel 2, Paganini et al. 2012; 
9P/Tempel 1, Mumma et al. 2005; 6P/d’ Arrest, Dello Russo et al. 2009). In short, it is critically 
important to develop a chemical taxonomy based on primary volatile composition. 

Infrared and radio spectroscopic observations of comets enable the measurement of gas 
production rates of primary volatile species. These surveys now permit quantitative sampling of 
eight or more primary species at infrared wavelengths, simultaneously with water, the (usually) 
dominant species. Several important primary volatiles (e.g., CH4, C2H2, C2H6, etc.) are readily 
measured at infrared wavelengths, but (being symmetric) their pure rotational transitions are 
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forbidden and so they cannot be observed at mm and sub-mm wavelengths. Radio observations 
are, however, very useful in studying the rotational transitions of more complex molecules 
(provided that they have a permanent dipole moment), and are thus complementary to infrared 
studies. 

Our group has identified ’’organics-enriched”, ’’organics-normal”, and ’’organics-depleted” 
comets on the basis of mixing ratios of organic volatiles, such as CH 3 OH, HCN, H 2 CO, C 2 H 2 , 
C 2 H 6 , CH 4 , and CO (Mumma et al. 2003; Crovisier et al. 2007; DiSanti and Mumma 2008). The 
Oort cloud comets C/2001 A2 (LINEAR) (Magee-Sauer et al. 2008), and C/2007 Wl (Boattini) 
(Villanueva et al. 2011a) have been identified as enriched in most organic species, and the Oort 
cloud comet C/1999 S4 (LINEAR) (Mumma et al. 2001; 2003) has been identified as severely 
depleted in most organics. Among ecliptic comets, 17P/HoImes is enriched (Salyk et al. 2007, 
Dello Russo et al. 2008) and 73P/Schwassman-Wachmann-3 is severely depleted (Villanueva et 
al. 2006, Dello Russo et al. 2007). However, 17P/Holmes was observed while outbursting at ~2.4 
AU (most comets are observed much closer to the Sun). It is likely that a fraction of the water in 
its inner coma was ejected in the form of ice (Yang et al. 2009), which (if depleted in organics) 
could cause enrichment in the organic gases measured at infrared wavelengths. Such polar and 
apolar ices were identified recently in comets 103 P/Hartley 2, C/2006 Wl (Boattini), and 
lOP/Tempel 2 (Mumma et al. 2011, Villanueva et al. 2011a, Paganini et al. 2012). Chemical 
diversity is found in each dynamical reservoir, and a more expansive sample of comets from both 
the Kuiper belt and the Oort cloud is needed in order to identify the true origin of an individual 
comet and the fractional representation of taxonomic classes within each dynamical reservoir. 

Here we present the organic composition of the ecliptic comet 2P/Encke, and discuss its organic 
depletion, its low rotational temperatures, and aspects of its possible origins. 2P/Encke (hereafter 
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Encke) is the prototypical “Encke-type” comet, with Tj = 3.025, period of 3.3 years, and 
aphelion distance of 4.09 AU. Encke‘s aphelion distance is unusually small, and dynamical 
simulations suggest that the comet may have been in a dormant state for a long period of time, 
after it became dynamically decoupled from Jupiter (Levison et al. 2006). Encke was discovered 
in 1786 by Mechain, at the Paris Observatory, and up to the time of our observations, it had 
returned to perihelion 67 times since its discovery (Sekanina 1991). 

2. OBSERVATIONS AND DATA REDUCTION 

Comet Encke reached closest approach to Earth (with A = 0.261 AU) on 17 Nov. 2003, and 
perihelion (q = 0.338 AU) on 30 Dec. 2003. We characterized the comet on 4-6 Nov. 2003 
(Table 1), using the cross-dispersed Near Infrared Echelle Spectrograph (NIRSPEC) at Keck-II 
atop Mauna Kea, HI. NIRSPEC has a 1024 x 1024 InSb detector array; and provides resolving 
power A/ AA ~ 25000 when using the 0.432" entrance slit. A single spectral frame (covering 1024 
x 1024 pixels) samples multiple echelle orders, e.g., orders 21-26 in the KL2 grating setting; and 
orders 22-27 in the KL1 grating setting (Table 1). This large spectral grasp is one of the main 
strengths of NIRSPEC since the simultaneous detection of multiple organic species and water 
minimizes numerous systematic uncertainties. The slit has a length of 24”, corresponding to 121 
pixels of 0.198” each in the spatial direction. 

We nodded the telescope by 12" (± 6" along the instrument slit) in a sequence Al, Bl, B2, A2 
(where A and B correspond to the two different positions of the comet in the slit). The difference 
in frames (A1-B1-B2+A2) and the further combination of the A and B beams cancels dark 
current and telescope thermal background, as well as sky emission lines (to second order in a 
Taylor series expansion about the mean air mass). We cropped individual echelle orders from the 
difference frame and (after flat field division) "masked" them in order to remove hot pixels and 
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cosmic ray hits. We re-sampled the data in the spatial and spectral directions, correcting for the 
initial tilt, to ensure that (in the re-sampled data) the pixels along each column correspond to a 
common frequency and that all pixels along a given row sample a common spatial position in the 
coma. 

We generated terrestrial transmittance spectra using our adapted (Villanueva et al. 2011b) 
version of the LBLRTM (Layer-by-Layer Radiative Transfer Model, Clough et al. 2005) spectral 
synthesis program, and convolved the fully resolved synthetic spectra to the instrumental 
resolution in order to determine the telluric abundances of H 2 O, CO 2 , O 3 , N 2 O, CO, C 2 H 6 and 
CH4 (by matching the synthetic and measured spectrum, through an iterative fitting process). We 
performed frequency calibration of the data by aligning modeled sky emission lines with those 
observed in the data frames. We scaled the convolved (to A./AA. ~ 25000) telluric model to the 
cometary continuum level; subtracting this isolated the residual emission spectrum of the comet, 
still multiplied by the atmospheric transmittance function. The fully-resolved best-fit model 
provided a precise value for the transmittance at each Doppler-shifted line position. We used 
spectra of a flux standard star to convert the counts per second in the comet data to flux density 
in W nT [cm' ]' . Detailed descriptions of our standard data analysis procedures are given in the 
literature (e.g., Radeva 2010, Villanueva et al. 2009, Bonev 2005, DiSanti et al. 2001, and 
references therein). 

3. RESULTS 

We present quantitative results for H2O, CH3OH, HCN, H2CO, C2H2, C2H6, CH4, and CO (Table 
2 ) in comet Encke. We present production rates for five species (H2O, CH3OH, HCN, C2H6, and 
CH4), 3 -a upper limits for three others (CO, H2CO, C2H2), and rotational temperatures for H2O, 


8 


CH 3 OH and HCN. We give mixing ratios for all trace species (production rates of organic 
species as percentages relative to water), for the three dates of observation. Results for Encke are 
compared with those for other comets in Table 3. 

3.1. Spectral Gallery 

Figure 1 presents a nucleus-centered, flux-calibrated spectral extract, which contains the spatially 
summed signal from nine spatial pixels centered on the nucleus, spanning a range of 1 .78" (~ 387 
km). The synthetic model of the terrestrial atmosphere is shown as a dashed line (superimposed 
on the cometary spectrum) in the upper panel of the figure. The corresponding residual spectrum 
of H 2 O (after subtraction of the transmittance model from the measured cometary spectrum) is 
shown in the lower panel, with the ±1 -ct noise envelope superimposed (dashed line). The noise 
envelope is largest at frequencies that correspond to sky emission lines (compare the upper and 
lower panels of Fig. 1). Figures 1-18 in the Electronic Supplementary Material present all 
remaining spectral orders for the three dates of observation. 

3.2. Measured Rotational Temperatures 


At infrared wavelengths, we sample the innermost coma over a range in nucleus-centered 
distance that depends on the comet’s geocentric distance (see Table 1). In this inner region, 
collisions thermalize the rotational populations (P m ) of the ground vibrational level of each 
molecule (Xie and Mumma 1992), and a rotational temperature (T rot ) is extracted to model these 


populations: Pm = 


g m exp(-Em/kTrot) 

Z(Trot) 


, where g m is the statistical weight of level 'm', E m is the 


energy of the level, and Z is the partition function. For Encke at the time of observation, the 
1.78” nucleus-centered extract spanned 387 km at the comet, corresponding to ± 194 km from 
the nucleus, as projected onto the sky plane. The number density of molecules at ~ 194 km from 
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the nucleus was approximately 8x10 cm' (for Q gas ~ 3 x 10‘ s' , see Discussion, §4). For a 
more active comet, the same number density would be measured at a greater distance from the 
nucleus, e.g., at ~ 346 km, even for the modestly productive comet C/2007 W 1 (Boattini) (Q ~ 
1.2 x 10 28 s' 1 ). Under these conditions, the rotational population of H 2 O in C/2007 W1 (Boattini) 
was well thermalized (retrieved T rot ~ 80 K) (Villanueva et al. 201 la), and in comet Encke the 
rotational population should be thennalized for all higher densities (i.e., within the sampled 
1.78” nucleus-centered extract). 

Accurate rotational temperatures are vital for calculating production rates of detected species, 
especially given that most measurements sample only a sub-set of ro-vibrational levels. To 
extract a rotational temperature for a given molecule, we compare the ratios of measured 
intensity and modeled fluorescence efficiency factor (g-factor) for the sampled spectral lines. 
The g-factors (in photons molecule ' 1 s' 1 ) predict the intensity of each line at a given T rot . 
Therefore, the temperature that most accurately describes the rotational level populations should 
provide the best agreement, within error, for the ratios of flux/g-factor for the measured lines. 
For the present analyses, we used the improved solar model and improved fluorescence models 
developed by our group (see Villanueva et al. 2011a (Boattini - solar & fluorescence models); 
2012a (H 2 O); 2012b (CH 3 OH), and original sources cited therein). 

The excitation analyses for three molecular species are shown in Fig. 2. They demonstrate our 
approach and the agreement (within error) among production rates obtained from individual lines 
at the retrieved (optimized) T rot , that corresponds to the solution with zero slope. The rotational 
temperatures derived for water, methanol, and hydrogen cyanide in Encke agree within 
experimental error (Table 2): T rot (H 2 O) = 29 ' 7 /+io K (4 Nov., K.L 1 setting), T rot (H 2 O) = 30 ' 2 /+ 2 
K (5 Nov., KL 1 setting), T r0 , (CH 3 OH) = 24 ' 5 /+6 K (5 Nov., KL1 setting), and T r0 , (HCN) = 20 ' 


10 


8 /+i6 K (6 Nov., KL2 setting). These temperatures are much lower than those we normally obtain 
for the inner coma of comets (see Discussion, §4). Fig. 3 illustrates the fit of the synthetic 
fluorescence models (convolved with modeled terrestrial transmittance) to the cometary spectra 
at each T rot . The rotational temperatures retrieved from FEO on the respective dates of 
observation were adopted for other species, for which T rot could not be extracted due to an 
insufficient spread of excitation energies amongst the measured lines. 

Our approach assumes that the observed medium is optically thin at infrared wavelengths. 
Optical depth effects must be considered for very active comets such as Ikeya-Zhang or Hale- 
Bopp (Dello Russo et al. 2004, DiSanti et al. 2001). In Encke, optical depth effects are negligible 
because the comet has very low gas production rates - (XFEO) is about 100 times smaller in 
Encke than in Ikeya-Zhang, and about 1000 times less than in Hale-Bopp. 

3.3 Production Rates 


We extracted a production rate Q, from each measured spectral line, assuming spherical outflow 


with uniform velocity. Then, Qi — ' 


4tcA Fi 


', where Fj is the flux of the i th line collected 


gil(hcu)tif(x) 

from a box on the sky centered on the nucleus and measuring 0.432” (slit-width) x 1.782” 
(extract length). Other parameters include the geocentric distance A [meters], the terrestrial 
transmittance (tj) at the frequency of the i th line; the fraction (f(x)) of the total coma content of 
the targeted species sampled by the beam; the energy (hcv) of a photon with wavenumber v (cm' 
') the line g-factor (gj) at T rot ; the molecular lifetime (t); and the assumed expansion velocity (0.8 
x Rh’ 0 5 km s' 1 ) (cf., Mumma et al. 2003). The weighted-mean production rate for each species is 
calculated from the individual production rates obtained from the sampled lines (using the 
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retrieved rotational temperature), where the weight is the inverse of the respective squared 
stochastic error for a given line (cf. Fig. 2). 

We first derive production rates from nucleus-centered extracts (which provide the highest 
signal-to-noise ratio), but the ‘nucleus-centered’ Q underestimates the ‘global’, or ‘terminal’, 
production rate owing to loss of flux in the instrumental slit (as a result of seeing effects or comet 
drift in the slit, etc.). Therefore, we determine the mean values of production rates (symmetric 
Qs) from extracts equidistant on both sides of the measured spatial intensity profile (Figure 4) - 
these reach a terminal value at ~ 1" from the nucleus. We calculate a growth factor for each 
species, as the ratio between the terminal and nucleus-centered production rate. Some molecules 
do not have sufficient signal-to-noise ratio away from the nucleus, which requires adopting a 
growth factor (usually from H2O, measured simultaneously, in the same setting), for the 
calculation of their terminal production rates. This procedure (averaging of production rates) also 
corrects for 1 -dimensional asymmetries in the gas outflow (Xie and Mumma 1996a, 1996b). 

We retrieved production rates for H 2 0.in the KL1, KL2, and MW_A setting, and for the 
following organic species: CH3OH, HCN, H 2 CO, C 2 H 2 , C 2 H6 (v 5 and v 7 bands), CH 4 and CO. 
We present 3-o upper limits for those species for which we cannot claim a firm detection: H 2 CO 
on 4 & 6 Nov.; C 2 F1 2 on 4 & 6 Nov.; C 2 H6 (V5 band) on 6 Nov., and CO on 5 & 6 Nov. CO 
measurements were combined from 5 & 6 Nov. to improve the signal-to-noise ratio, given the 
limited time on source in the MW_A setting. The instrumental setting, time on source, global 
production rate and mixing ratio for each molecule is presented in Table 2. The production rates 
of the dominant volatile - water, agree within 2 -ct on all three dates, and for all three settings 
used. 
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DISCUSSION 


Our high-resolution infrared spectra of comet Encke on three consecutive dates reveal its unusual 
composition, and very low rotational temperatures consistent with its low water production rate. 
6P/d’ Arrest is another example of a comet with a very low production rate (comparable to that of 
Encke), which is also characterized by low rotational temperatures (30-50 K) (Dello Russo et al. 
2009). The local collisional excitation rate (c) is directly proportional to the production rate of a 

comet, according to: c = nav t , where n is the number density (n ^—7), a is the collisional 

47tv 0 r" 

cross-section for excitation, v t is the thermal velocity, and v 0 is the gas outflow velocity (Weaver 
and Mumma 1984). Is it possible that the low rotational temperatures observed in comet Encke 
are the signature of relaxation by radiative decay dominating collisions in controlling its 
rotational level population? 

Xie and Mumma (1992) examined this question in a detailed study of water in comet I P/Halley 
that included collisions of electrons with water, water with water, infrared fluorescent pumping 
of H2O with rotational cascade, and radiative decay of H2O rotational levels. They applied their 
model for conditions measured during the Giotto coma flythrough (14 March 1986; 0.9 AU). 
For a water production rate of 5 x 10“ s' and kinetic temperature of 60 K, they evaluated the 
excitation & de-excitation rates over the nucleocentric distance range of 1000 km to 3 x 10' km. 
In their model, the water density was ~ 6 x 10 6 cm° at 3000 km, well within the coma region 
where FEO rotational emission lines are optically-thick and radiative trapping reduces the 
escape-to-space probability significantly. Xie and Mumma (1992) showed that electron and 
neutral-neutral collisions dominate rotational excitation and de-excitation of water under such 
conditions. Since cross-sections and radiative decay times are similar for other polar gases 
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(CH 3 OH, H 2 CO, HCN, CO), we expect they too will be controlled by collisions. Radiative decay 
rates are zero for apolar gases (CH 4 , C 2 H 6 , etc.), so collisions will control their rotational 
temperatures so long as fluorescent pumping rates remain small. This is true within the volume 
sampled in our investigation, since the pumping rates for all molecules are less than one event 
per 1000 seconds. For conditions in 2P/Encke (Q = 2.6 x 10 27 s' 1 , v 0 ~ 0.8 km/s), the density of 6 
x 10 6 cm" ' is reached at ~200 km nucleocentric distance. Most of the molecules sampled in our 
beam fall within this distance range. We therefore conclude that rotational populations in all 
molecules sampled in this study (Table 2) are controlled by collisions. 

In an attempt to explain the high rotational temperatures found in the near nucleus region of 
73P/Schwassmann-Wachmann 3, Fougere et al. 2012 suggest that sublimation of water from 
warm grains introduced excess heat to the near-nucleus coma, thus preventing adiabatic cooling 
of the expanding gas and maintaining rotational temperatures far above the adiabatic limit. A 
similar effect is consistent with the release of water-rich icy grains in C/2006 W 1 (Boattini), as 
inferred from the distinct spatial distributions of polar and apolar volatiles in the coma 
(Villanueva et al. 2011a). Moreover, our observed (unusual) depletions of primary volatiles in 
Encke, and its unusual thermal history, suggest that release of ice-mantled grains (and the 
concomitant excess heat injected by sublimation) is not likely. Thus, we may be seeing the 
expected adiabatic cooling of expanding gas. 

We must also consider the effects of low thermalization efficiency of fast H-atoms, which are 
by-products of H 2 O and OH photodissociation, and whose thermalization increases the 
temperature of the coma (for comets with high gas production rates) after the initial temperature 
decrease introduced by expansive adiabatic cooling (Combi et al. 2004). However, re-heating 
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normally operates at greater nucleocentric distances than we sampled in Encke. A detailed 
evaluation of these competing effects is deferred to a future publication. 

Historically, Encke’s spectral continuum at optical wavelengths has long been noted as unusual. 
Optical observations demonstrate a very weak continuum consistent with a deficiency of sub- 
micrometer sized dust in Encke. This behavior and observations of a strong thermal IR 
continuum (Ney 1974, Campins et al. 1982) were argued to imply a lower than typical fraction of 
small dust particles in Encke’s coma (i.e., most dust particles had radii between 5 and 10 pm, or 
even > 20 pm) (Gehrz & Ney 1989, Jewitt et al. 2004, Lisse et al. 2004). We observed a strong 
infrared continuum, consistent with a substantial contribution from micron-sized grains. 

The gaseous activity of Encke reported here agrees with contemporaneous measurements 
obtained from the Odin space observatory. Our water production rate in Encke of (2.6 ± 0.1) x 
10 s" (weighted mean from all settings, 4-6 Nov.) at a mean heliocentric distance of 1 .2 AU, is, 
as expected, lower than that measured closer to the Sun: Q(H 20 ) = (4.9 ± 0.7) x 10 27 s' 1 , at Rh of 
1 .01 AU (16 Nov. 2003), from Odin observations of the 557 GHz line of H 2 O (Biver et al. 2007) 
(although there is agreement at the 3-a level). Scaling our production rate, to account for 
differences in heliocentric distance, by Rh" (Rh = 1.2 AU), provides agreement at the 1-a level 
with Odin measurements. 

Encke is an important addition to the growing database of ecliptic comets that have been 
characterized at infrared wavelengths. Its composition is unusual in being depleted in most 
organic species (HCN, H 2 CO, C 2 H 2 , C 2 H 6 , CH 4 , and CO) but enriched in CH 3 OH. The weighted- 
mean mixing ratios for measured species are presented in Table 3, along with those of other 
Jupiter-family comets (21P/Giacobini-Zinner, 73P/Schwassman-Wachmann-3, 9P/Tempel 1, 
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6 P/d’ Arrest, lOP/Tempel 2 and 17P/Holmes) and three Oort cloud (current) ‘end-member’ 
comets (the severely depleted C/1999 S4 (LrNEAR), and the significantly enriched comets 
C/2001 A2 (LINEAR) and C/2007 W1 (Boattini)). The only organics-enriched member of the 
presented Jupiter-family group is 17P/Holmes, whose enrichment is most likely overestimated, 
given its heliocentric distance of observation (2.4 AU) - where water in the inner coma is more 
likely in the form of ice (as previously discussed (§1)). Table 3 indicates that: (1) the number of 
ecliptic comets observed at infrared wavelengths is insufficient for statistical analysis, and no 
clear trends are seen in the current sample; (2) Encke is the ecliptic comet with highest methanol 
mixing ratio, studied in the infrared. 

Encke’s composition is also presented in Fig. 5, in comparison to that of the severely depleted 
JFC 73P/Schwassman-Wachmann-3C (hereafter 73P/S-W 3-C), the severely depleted C/1999 S4 
(LINEAR) (hereafter C/1999 S4), and the enriched C/2001 A2 (LINEAR) (hereafter C/2001 A2) 
and C/2007 W1 (Boattini). Fig. 5 clearly illustrates the significant depletion of Encke in HCN, 
H 2 CO, C 2 H 2 & CO; its moderate depletion in C 2 H 6 and CH 4 ; and its enrichment in CH 3 OH. 

Our weighted mean mixing ratios (from different dates) for CH 3 OH (3.48 ± 0.27%), H 2 CO (< 
0.13% at 3 -ct), and HCN (0.09 ± 0.01%) in Encke are consistent with those measured at radio 
wavelengths: CH 3 OH (4.1%), H 2 CO (< 1.4 %, 3-ct), and HCN (0.09%) (Crovisier et al. 2007). 
Infrared and radio measurements both indicate enrichment of Encke in CH 3 OH and depletion in 
HCN, and our infrared measurements alone indicate depletion in CH 4 , C 2 H 2 , C 2 H 6 , H 2 CO and 
CO. By contrast, optical observations of a sample of 85 comets found Encke to be ‘typical’ in 
carbon-chain molecules (A’Hearn et al. 1995; see also Fink et al. 2009), but Cochran et al. 
(2012) find Encke to be CH/CN-enriched. However, the product species sampled at optical 
wavelengths can have multiple origins (including release from grains), which greatly complicates 
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drawing cosmogonic implications from their measured abundance ratios. One possible 
explanation for the discrepancy between infrared (confirmed by radio) and optical measurements 
is that some carbon-chain parent molecules not sampled in the infrared may be more abundant in 
Encke, resulting in typical values for C 2 and C 3 but depletion in C 2 H 6 and C 2 H 2 (acetylene is a 
possible parent of C 2 , but C 2 H 6 is not (Weiler 2012)). Several other comets also show discordant 
comparisons in product and primary volatile (i.e., potential parent) species (§ 1 ). 

It is very difficult to determine whether the depletion of 2P/Encke is evolutionary or primordial. 
Successive depletion during its numerous orbits around the Sun could have induced preferential 
loss of the most volatile ices (CO and CH 4 , having sublimation temperatures of 25 K and 31 K, 
respectively) over those having higher sublimation temperatures (e.g., CH 3 OH, 99 K) 
(Yamamoto 1985). For example, ecliptic comet 73P/S-W 3-C, has a short period (5.4 years) and 
is a severely depleted end-member. However, the ecliptic comet 103P/Hartley 2 (P = 6.5 years) 
is normal in its organic primary composition (although depleted in C 2 H 2 , and highly enriched in 
C0 2 (A’Heam et al. 2011)), and ‘typical’ in its product species. Taken together, these 3 
examples suggest that something other than evolutionary processing (e.g., formative conditions) 
influenced current volatile compositions, although Encke with its very short orbital period and 
small perihelion distance may represent a special case. 

It is possible that if an ecliptic comet formed later in the Solar System’s evolutionary timeline 
than did Encke, more ionizing x-ray flux from the Sun would reach greater distances thanks to 
nebular clearing, and the resulting higher H-atom densities - combined with lower temperatures 
enhancing H-retention on the surfaces of grains - would increase the efficiency of H-atom 
addition reactions (e.g., more C 2 H 6 would form from C 2 H 2 , and CH 3 OH from CO; see below). 
The conversion efficiency for C 2 H 2 into C 2 H 6 can be calculated as + C 2 H 6 ), and is > 
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0.8 in Encke. This relatively high value could result from formation in a region with higher El- 
atom densities, and/or lower temperatures, than were experienced by comets with a lower CiEE 
hydrogenation efficiency. High hydrogenation efficiency is also consistent with the low upper 
limits for CO & H 2 CO, but significant enrichment of CH 3 OH in Encke: (H 2 CO + CH 30 H)/(C 0 
+ H 2 CO + CH 3 OH) > 0.7. Alternatively, a high initial endowment of CO with low conversion to 
H 2 CO and CH 3 OH could indicate that the 'missing’ CO (the most volatile species) was lost 
during Encke’s numerous orbits around the Sun. The relatively efficient conversion implied for 
C 2 H 2 (to ethane) may also suggest that most CO that condensed onto its pre-cometary grains 
later converted to CH 3 OH. If true, the severe depletion of CH 3 OH in C/1999 S4 together with its 
enrichment in Encke is consistent with this view. 

It is important to note that comet Encke is not as severely depleted as C/1999 S4, which comes 
from the Oort cloud and is considered to be dynamically new. It was previously suggested that 
C/1999 S4 is depleted because it likely formed closer to the young Sun than ‘organics-normal’ 
comets (Mumma et al. 2001). That is a possible explanation for comet Encke’s depletion as well, 
but if so it probably formed farther from the Sun than did C/1999 S4, given Encke’s less severe 
depletion. 

In addition to heliocentric distance, other factors must be considered to explain the depletion (or 
enrichment) of a cometary nucleus, such as localized heating, or formation of the nucleus from 
fractions originating in regions with diverse chemical composition. Our quantitative production 
rates on three consecutive dates can test for potential heterogeneity of Encke’s nucleus. If the 
nucleus were heterogeneous, we would expect to see variability in the measured mixing ratios as 
it rotates and different vents are sequentially activated under the influence of solar radiation 
(depending on the size of the fractions, locations of vents, etc.). The rotation period of Encke was 
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measured to be 11.056 ± 0.024 h (sidereal, Jockers et al. 2011), 11.083 ± 0.003 h (synodic, 
Lowry and Weissman 2007), and 1 1.079 ± 0.009 h (synodic, Fernandez et al. 2005). Our KL2 
observations were separated by ~ 47 h 30 m (4.3 rotation periods) and we likely observed release 
from different regions of Encke’s nucleus, but our KL1 observations were separated by ~ 22 h 
(2.0 periods), and may sample the same region. We present a comparison of the mixing ratios in 
comet Encke on different dates in Figure 6, and see agreement in the abundances measured for 
each species within 1-a (2-a for F1CN). This lack of temporal variability in the measured organic 
abundances is consistent with compositional homogeneity in the mean volatile release (averaged 
over all active vents, regardless of possible chemical diversity), similar to comet 103P/Hartley-2 
(Mumma et al. 2011, Dello Russo et al. 2011), C/2000 WMi (Radeva et al. 2010), and C/2004 
Q2 (Bonev et al. 2009). However, our sparse temporal sampling and limited spatial information 
prevent us from drawing a definitive conclusion on this question. 

5. SUMMARY 

The ecliptic comet Encke is an important addition to our growing database of comets studied at 
infrared wavelengths, first, because fewer ecliptic comets have been studied in comparison to 
Oort cloud comets, and second, because it is a dynamical end-member owing to its very short 
orbital period of 3.3 years. Our investigation of Encke resulted in the following conclusions: 

1. Encke exhibits one of the lowest gas production rates for water among comets in our 
sample, with Q(H20) being (2.6 ± 0.1) x 10 27 s" 1 (weighted mean from all settings, 4-6 
Nov.) at a mean heliocentric distance of 1.2 AU. The scatter among individual 
measurements is somewhat larger than the error in the mean (Table 2). 
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2. Encke exhibits very low rotational temperatures in the near-nucleus coma (T ro t ~ 20 - 30 
K, for H2O, CH3OH and HCN). Collisions are sufficient to maintain the rotational level 
populations. The absence of excess heating by water released from icy-mantled grains 
likely permits the predicted adiabatic cooling of outflowing gas in the near nucleus coma 
to control the rotational temperature. 

3. Encke is depleted in HCN, H2CO, C2H2, C2H6, CH4 and CO, but is enriched in CH3OH 
(it being the least volatile primary species other than H2O). Our measurements of 
CH3OH, H2CO and HCN agree with those from radio observations. The comet may have 
lost most of its volatiles during its numerous orbits around the Sun, or it may have 
formed closer to the young Sun than did organics-normal and organics-enriched comets. 

4. Mixing ratios in Encke are consistent on different dates, but the temporal sampling of our 
data is insufficient for a firm conclusion as to whether its nucleus is homogeneous in the 
bulk. More observations of Encke are necessary to confirm its homogeneity. 

In our growing database of Oort cloud and ecliptic comets, we see no correlation of chemical 
composition with originating dynamical reservoir or orbital period. For example, the 
dynamically new comet C/1999 S4 is more severely depleted than Encke, and so is the 5.4 year 
period (split) ecliptic comet 73P/SW 3-B and -C; while the 6.5 year period ecliptic comet 
103P/Hartley 2 is normal in its composition (except for its depletion in acetylene). More comets 
need to be observed at infrared wavelengths, which sample primary volatiles native to the 
cometary nucleus, and uniquely sample the symmetric hydrocarbons. The ‘organics-depleted’, 
‘organics-normal’ and ‘organics-enriched’ (along with ‘atypical’) compositional classes should 
be regarded as a (plausible) working hypothesis. A statistically significant sample of comets will 
permit us to test this classification. At present we do not know if comets currently categorized as 
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“depleted” or “enriched” are in fact compositional end members. Future observations would 
elucidate whether the current classes are well separated or if there is a more continuous 
distribution of abundances amongst the comet population. Measurements of cosmogonic 
parameters, such as spin temperature of water (indicating its formation temperature) and D/H 
ratios in comets, are needed in order to evaluate cometary composition as one of the keys in 
understanding the formation and evolution of our Solar System. 
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Table 1: Log of Observations of comet Encke. 


UT Time (2003) 

Integration 
time [min] 

Setting 1 

Frequency range 
of setting [cm 1 ] 

Rh 11 

[AU] 

A" 

[AU] 

A 11 
A dot 

[km s" 1 ] 

4 Nov., 5:09-7:17 

84 

KL2 

3464-2759 

1.210 

0.313 

-13.59 

9 : 29 - 9:52 

20 

KL 1 

3521 -2833 

1.208 

0.312 

- 13.19 

5 Nov., 7:23-8:30 

48 

KL1 

3521 -2833 

1.194 

0.305 

-12.45 

9:31-10:05 

12 

MW A 

2162- 1996 

1.193 

0.304 

-12.25 

6 Nov., 4:40-6:21 

80 

KL2 

3464-2759 

1.180 

0.299 

-11.78 

9:18-9:50 

14 

MW A 

2162- 1996 

1.178 

0.297 

-11.31 


‘The KL1 setting simultaneously samples CH 3 OH, C 2 H 6 & H 2 O. The KL2 setting 

simultaneously samples HCN, CH 4 , C 2 H 2 , C 2 H 6 , H 2 CO & H 2 O; and the MW_A setting 
simultaneously samples CO & H 2 O. 

M Rh is the heliocentric distance, A is the geocentric distance, and Ad 0 t is the line-of-sight 
velocity. 
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Table 2. Production rates and mixing ratios of volatiles in 2P/Encke.' 


Setting / Time on 

Molecule (T rot ) 

Global Q 

Mixing Ratio % 

Source 


[10 25 s' 1 ] 



4 Nov. 2003 


KL 1 / 20 min 

H 2 0 

(29 - 7 / +10 K) 

215.7 ± 25.3 

100.00 

C 2 H 6 V7 
(29 K) 

0.7 ±0.1 

0.32 ±0.05 

CH 3 OH 
(29 K) 

8.5 ± 1.5 

3.93 ±0.69 

KL2 / 84 min 

H 2 0 
(29 K) 

272.2 ± 86.7 

100.00 

H 2 CO 
(29 K) 

< 0.4 at 3a 

< 0.13 at 3a 

C 2 H 6 V5 
(29 K) 

1.0 ±0.4 

0.36 ±0.09 

CH 4 
(29 K) 

0.9 ±0.4 

0.31 ±0.11 

HCN 
(29 K) 

0.4 ±0.1 

0.13 ±0.03 

c 2 h 2 

(29 K) 

< 0.3 at 3a 

<0.10 at 3a 

5 N 

ov. 2003 

KL1 / 48 min 

H 2 0 

(30 ' 2 / +2 K) 

298.0 ±23.9 

100.00 

C 2 H 6 V7 
(29 K) 

0.9 ±0.1 

0.31 ±0.04 

CH 3 OH 
(24 - s / +6 K) 

9.1 ± 1.0 

3.01 ±0.28 

CH 3 OH 
(30 K) 

10.2 ± 1.0 

3.40 ±0.29 


MWA Combined 

12 min (Nov. 5) 
& 

14 min (Nov. 6 )' 

H 2 0 
(30 K)“ 

321.8 + 82.0 

100.00 

CO 

(30 K ) 11 

< 5.7 at 3a 

< 1.77 at 3a 

6 IN 

ov. 2003 

KL2 / 80 min 

H 2 0 
(30 K) 

256.4 + 20.4 

100.00 

H 2 CO 
(30 K) 

< 0.3 at 3a 

< 0.13 at 3a 

C 2 H 6 V 5 
(30 K) 

< 0.6 at 3a 

< 0.23 at 3a 

ch 4 

(30 K) 

1.3 ±0.7 

0.52 + 0.28 

HCN 

(20 - 8 / +16 K) 

0.18 ± 0.03 

0.07 + 0.01 

HCN 
(30 K) 

0.21 ±0.04 

0.08 + 0.01 

C 2 H 2 
(30 K) 

< 0.2 at 3a 

< 0.08 at 3a 


'Production rates for H 2 O assume an equilibrium Ortho-to-Para ratio of 3.0. Detailed 
analysis of the OPR in comet Encke is reserved for a further publication. 


"The production rates of H 2 O and CO for MW_A are based on data from 5 Nov. (12 
min on source) and 6 Nov. (14 min on source), which were combined to provide a 
higher signal to noise ratio. 


Table 3. Mixing ratios in 2P/Encke, in other Jupiter-family comets; and in ‘organics-depleted’ and 
‘organics-enriched’ end-member Oort Cloud comets. 


Mixing Ratios in a Sample of Jupiter-family Comets 


Comet 

CH 3 OH 

HCN 

H 2 CO 

c 2 h 2 

c 2 h 6 

ch 4 

CO 

73P/S-W 3-C 1 

0.149 
± 0.029 

0.242 

±0.014 

0.147 
± 0.033 

0.049 
± 0.020 

0.107 

±0.011 

<0.25 

0.5 

±0.13 

2 1 P/Giacobini-Zinner 11 

1.22 

±0.11 

<0.27 

<0.8 

<0.42 

0.136 
± 0.023 

- 

10.25 

±5.80 

9P/Tempel l 111 
(pre-impact) 

1.3 

±0.2 

0.18 

±0.06 

- 

- 

0.23 

±0.04 

- 

- 

9P/Tempel l 111 
(post-impact, 
including ejecta) 

0.99 

±0.17 

0.22 

±0.03 

- 

0.13 

±0.04 

0.40 

±0.04 

0.54 

±0.30 

4.3 
± 1.7 

6P/d’ Arrest lv 

1.42 

±0.30 

0.034 
± 0.008 

0.36 

±0.09 

<0.052 

0.256 

±0.061 

- 

- 

lOP/Tempel 2 V 

1.58 

±0.23 

0.13 

±0.02 

<0.11 

<0.07 

0.39 

±0.04 

- 

- 

103 P/Hartley 2 VI 

2.28 

±0.38 

0.26 

±0.02 

0.23 

±0.05 

0.08 

±0.01 

0.75 

±0.07 

< 1.2 


1 7P/Holmes v " 

3.0 
± 1.0 

0.538 
± 0.075 

- 

0.344 
± 0.053 

1.90 

±0.54 

- 

- 

2P/Encke vin 

3.48 

±0.27 

0.09 

±0.01 

<0.13 

<0.08 

0.32 

±0.03 

0.34 

±0.10 

< 1.77 


Mixing Ratios in End-Member Oort Cloud Comets 


C/1999 S4 IX 
(Organics- Depleted) 

<0.15 

0.10 

±0.03 

- 

<0.12 

0.11 

±0.02 

0.18 
± 0.06 

0.9 

±0.3 

C/2001 A2 X 

3.9 

0.6 

0.24 

0.5 

1.7 

1.2 

3.9 

(Organics-Enriched) 

±0.4 

±0.1 

±0.05 

±0.1 

±0.2 

±0.2 

± 1.1 
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C/2007 W1 (Boattini) Xi 

3.69 

0.50 

< 0.12 

0.29 

1.97 

1.57 

4.50 

(Organics-Enriched) 

±0.07 

± 0.01 


± 0.02 

±0.03 

±0.16 

±0.51 


'Dello Russo et al. 2007 (C 2 H 6 , C 2 H 2 , HCN, H 2 CO & CH 3 OH, 14.5 May 2006), Villanueva et al. 
2006 (CH 4 , 7 Apr. 2006), DiSanti et al. 2007a (CO, 27 & 30 May 2006). 

"Weaver et al. 1999 (C 2 H 2 , HCN, H 2 CO, Oct. 25-29 1998), Mumma et al. 2000 (CO, Oct. 2-10 
1998) & DiSanti et al. 2012 (C 2 H 6 & CH 3 OH, Jun. 3 2005). 

'"Mumma et al. 2005 (Jun. 3 2005 & Jul. 4-5 2005), DiSanti et al. 2007b (4 Jul. 2005), DiSanti and 
Mumma 2008. 

IV Dello Russo et al. 2009 (1 1 Aug. 2008). 
v Paganini et al. 2012 (26 Jul. 2010). 

vl Mumma et al. 201 1 (22 Oct. 2010, mixing ratios in agreement with those presented for other 
dates; CH 4 is for 16 Nov.). 

vll DeIlo Russo et al. 2008 (27 & 3 1 Oct. 2007). 

vlll This work. Values are weighted means from different dates; 3-c upper limits are provided where 
indicated. 

IX Mumma et al. 2003 (13 Jul. 2000). 

x Magee-Sauer et al. 2008 & Gibbet al. 2007 (9 & 10 Jul., 4 & 10 Aug. 2001). 

XI 

Villanueva et al. 201 la (9 & 10 Jul. (weighted means), mixing ratios relative to symmetric water). 
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Wavenumber |cm '| 


Figure 1: 2P/Encke cometary spectrum and superimposed terrestrial transmittance model 
(dashed line) for KL1, Order 27, on 5 Nov. 2003 (upper panel), and residual spectrum (lower 
panel). Spectral lines of H 2 O, OH, and blends are seen. The flux density scale is shown at left, 
and the (±l-o) noise envelope is shown as dashed lines centered on zero flux density, in each 
panel. 
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a 



Excitation Energy [cm 0 ] Excitation Energy [cm 1 ] Excitation Energy [cm 1 ] Excitation Energy [cnr 1 ] 


Figure 2: Excitation analysis for H2O on 4 Nov. (Panel A), CH3OH on 5 Nov. (Panel B), H2O 
on 5 Nov. (Panel C), and HCN on 6 Nov. (Panel D). At the retrieved T rot , the agreement among 
production rates derived from individual measured lines is satisfactory, for each species. 


Rotational temperatures for H2O in the KL1 setting (Panels A & C) are retrieved from the 
combination of spectral lines in echelle orders 26 and 27. 
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Figure 3: Residual features for HiO on 4 Nov. (Panel A 1 for Order 26 & Panel A2 
for Order 27); CH3OH on 5 Nov. (Panel B), H2O on 5 Nov. (Panel Cl for Order 
26 & Panel C2 for Order 27), and HCN on 6 Nov. (Panel D), with the fluorescence 
model for each molecule superimposed in red. The combination of spectral lines 
sampled in orders 26 and 27 for H2O on 5 Nov. resulted in the best-constrained 
T ro t of 30 ' 2 /+ 2 K. 
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Nucleus Nucleus 
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Figure 4: Emission spatial profiles for H 2 O and C 2 H 6 (KL1 setting, Nov. 5) in black, and dust 
continuum in red. Asymmetry in the spatial profile is compensated by extracting production rates 
on both sides of the profile, and reporting an averaged global production rate. The continuum 
profile is clearly displaced from the gas profile, suggesting preferential outgassing towards the 
sun. The orientation of the instrument slit (89°) is shown on the sky plane, as well as the aspect 
in the Sun-Nucleus-Observer plane. 
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Volatile Abundances Relative to H z O 



Figure 5: Comparison among mixing ratios of organic species in the severely depleted 
ecliptic comet 73P/S-W 3 -C, the severely depleted Oort cloud comet C/1999 S4, the 
ecliptic comet Encke, and the enriched Oort cloud comets C/2001 A2 and C/2007 W 1 
(Boattini). Organic species are presented on the figure in terms of decreasing sublimation 
temperature (from left to right). Comet Encke is enriched in CH 3 OH; and depleted in 
HCN, H 2 CO, C 2 H 2 , C 2 H 6 , CH 4 & CO. 


Mixing Ratio [%] 


Volatile Abundances Relative to H 2 0 



Figure 6: Agreement among mixing ratios of organic species in 2P/Encke on 4-6 Nov. 
2003, within 1-a (2-c for HCN). The lack of variability would suggest homogeneity of 
Encke’s nucleus, however, we lack sufficient temporal sampling to make a firm 
conclusion. 
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